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We calculate the incoherent resonant and non-resonant scattering production of sterile neu- 
trinos in the early universe. We find ranges of sterile neutrino masses, vacuum mixing angles, 
and initial lepton numbers which allow these species to constitute viable hot, warm, and cold 
dark matter (HDM, WDM, CDM) candidates which meet observational constraints. The con- 
straints considered here include energy loss in core collapse supernovae, energy density limits 
at big bang nucleosynthesis, and those stemming from sterile neutrino decay: limits from ob- 
served cosmic microwave background anisotropies, diffuse extragalactic background radiation, 
and 6 Li/D overproduction. Our calculations explicitly include matter effects, both effective mix- 
ing angle suppression and enhancement (MSW resonance), as well as quantum damping. We 
for the first time properly include all finite temperature effects, dilution resulting from the an- 
nihilation or disappearance of relativistic degrees of freedom, and the scattering-rate-enhancing 
effects of particle-antiparticle pairs (muons, tauons, quarks) at high temperature in the early 
universe. 
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I. INTRODUCTION 

In this paper we calculate the non-equilibrium resonant 
and non-resonant production of sterile neutrinos in the 
early universe and describe cosmological and astrophysi- 
cal constraints on this dark matter candidate. Depending 
on their masses and energy distributions, the sterile neu- 
trinos so produced can be either cold, warm, or hot dark 
matter (CDM, WDM, HDM, respectively) and may help 
solve some contemporary problems in cosmic structure 
formation. We can define sterile neutrinos generically as 
spin-1/2, S'?7(2)-singlet particles which interact with the 
standard S'C/(2)-doublet ("active") neutrinos u e , v^, and 
v r , solely via ordinary mass terms. Singlet neutrinos with 
masses ~ 10 12 GeV arise naturally, for example, in "see- 
saw" models of neutrino mass in grand unified theories 
(GUTs) |Q. Recent solar, atmospheric, and accelerator 
neutrino oscillation experiments P, ||[ ^[ ||), however, 
imply the existence of four light neutrino species with 
masses < 10 eV, only three of which can be active, on 
account of limits on the invisible decay width of the Z° 
boson @. The remaining neutrino must be sterile. The 
existence of multiple generations of quarks and leptons 
in the standard model (SM) of particle physics, as well as 
many independently- motivated extensions of the SM, im- 
ply that there are additional, more massive sterile neutri- 
nos with couplings to active neutrinos currently beyond 
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direct experimental reach. We describe herein the cos- 
mological implications of these heavier sterile neutrinos. 

Some early constraints on massive sterile neutrino pro- 
duction in the early universe were derived in Refs. [g, 
The first analytical estimates of the relic sterile neu- 
trino abundance from scattering-induced conversion of 
active neutrinos were made by Dodelson and Widrow 
|l0| . They assumed a negligible primordial lepton num- 
ber, or asymmetry, so the neutrinos are produced non- 
resonantly. They found that sterile neutrinos could be a 
WDM candidate with interesting consequences for galac- 
tic and large scale structure. A process that can create 
sterile neutrino dark matter with a unique energy dis- 
tribution was proposed by Shi and Fuller O]. In that 
work a non-vanishing initial primordial lepton number 
gives rise to mass-level crossings which enhance sterile 
neutrino production, yielding dark matter with an en- 
ergy spectrum which is grossly non-thermal and skewed 
toward low energies. In particular, the average energy is 
significantly less than that of an active neutrino, and the 
neutrinos suppress structure formation on small scales 
and behave like CDM on large scales (i.e., they behave 
as WDM or "cool DM"). 

At no previous time has there been a greater need for 
a more comprehensive study of the cosmological and as- 
trophysical consequences of active-sterile neutrino mix- 
ing. The continuing influx of data from neutrino experi- 
ments and new observations of galaxy cores and clusters 
demand a more detailed understanding of the physics of 
sterile neutrino dark matter production and more sophis- 
ticated calculations of their relic abundances. 

We examine these experimental and observational is- 
sues in Sees. |H] and [II, respectively. In Sec. [V we discuss 
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the present limits on primordial lepton asymmetries and 
possible dynamical origins of these asymmetries. We at- 
tempt to unify the perspectives of previous work, so we 
take the initial lepton asymmetry of the universe to be a 
free parameter within the rather generous bounds set by 
experiment. After reviewing the physics of neutrino os- 
cillations and matter-affected neutrino transformation in 
Sec. [v| we compute the consequences of pre-existing lep- 
ton asymmetries f or st erile neutrino dark matter scenar- 
ios in Sees. VI and VII. In Sec. VI, we give the Boltzmann 
equations governing the non-equilibrium conversion of 
active neutrinos into sterile neutrinos and solve them in 
the limit of non-resonant conversion, a limit which ob- 
tains when the lepton asymmetry is su fficie ntly small, of 
order the baryon asymmetry. In Sec. VII, we consider 
larger asymmetries and examine in detail the physics of 
resonant neutrino transformation in the early universe. 
For the non-resonant and several resonant cases we dis- 
play contours of constant relic density of sterile neutri- 
nos in the plane of neutrino mixing parameters. These 
results confirm and extend previous work For 
both sets of calculations we also take into account finite- 
temperature and finite-density effects for all three active 
neutrino species; the dilution of sterile neutrino densities 
due to heating of the photons and active neutrinos from 
the annihilation of particle-antiparticle pairs before, dur- 
ing, and after the quark-hadron (QCD) transition, and 
we allow for the enhanced active neutrino scattering rates 
due to the increased number of sca tterers in equilibrium 
at high temperatures. In Sec. VIII we calculate the colli- 
sionless damping scales relevant for structure formation, 
scales which classify the regions of dark matter param- 
eter space as HDM, WDM, or CDM regions. We con- 
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from the diffuse cxtragalactic background radiation (DE- 
BRA), cosmic microwave background (CMB), big bang 
nucleosynthesis (BBN), and 6 Li and D photoproduction. 
Finally, in Sec. H we examine the implications of active- 
sterile neutrino mixing in core-collapse (Type Ib/c, II) 
supernovae. Conclusions are given in Sec. XI. Through- 
out the paper we use natural units with h = c = fcs = 1. 



II. NEUTRINO ANOMALIES 

Recent experiments have provided data indicating ev- 
idence for new neutrino physics. The most significant 
recent evidence is the Super-Kamiokande Collaboration's 
statistically convincing result 0] , verifying previous mea- 
surements 0, of a suppression of the atmospheric v^jv^ 
flux. The most persuasive Super-Kamiokande evidence 
for neutrino oscillations is the measured zenith angle de- 
pendence of the Vpjv^ flux, an observation fit most sim- 
ply by maximal v ^ ^ v T mixing in vacuum, with vacuum 
mass-squared difference Sm 2 ~ 3 x 10 _3 eV 2 . 

On another front, the ground-breaking observations by 
the Homestake Collaboration found a solar neutrino flux 
far below that predicted on the basis of sophisticated 



solar models j|. The solar neutrino problem has an in- 
teresting possible solution through matter-enhanced res- 
onant conversion via the Mikheyev-Smirnov-Wolfenstein 
(MSW) mechanism p2[ or through vacuum or "quasi- 
vacuum" oscillations ]13| . Depending on whether the so- 
lar solution involves two-, three-, or four-neutrino mixing, 
the parameter space of neutrino mass-squared difference 
and vacuum mixing angle are constrained differently 

A third indication for neutrino oscillations comes from 
the Los Alamos Liquid Scintillator Neutrino Detector 
(LSND) Collaboration's observations of excess v e and 
v e events in beams of and respectively. These 
have been interpreted as evidence for neutrino oscilla- 
tions in the — ► v e and — > v e channels ||. The 
Karlsruhe Rutherford Medium Energy Neutrino (KAR- 
MEN) experiment probes the same channels but does not 
see evidence for neutrino oscillations jl4|. A joint analy- 
sis of the LSND and K ARMEN data has found that there 
are regions of neutrino mixing parameter space consistent 
with both experiments' results 

Efforts have been made to embed the above neutrino 
oscillation solutions within a three-neutrino framework 
Jl6| , [It], [l^, ^9|. Leaving aside maximal vacuum mixing 
of all three neutrino species [[t6|, these analyses gener- 
ally require the atmospheric neutrino oscillation length 
scale to be the one associated with the short-base-line 
LSND experiment jlT], |l8| , or the solar and atmospheric 
solutions must be built on the same mass difference scale 
& 

However, the zenith-angle dependence of the Super- 
Kamiokande measurement requires the atmospheric neu- 
trino oscillation length to be much larger than the cor- 
responding LSND scale q|. This disfavors the first 
three-neutrino scheme. Additionally, the three solar 
neutrino experimental modes presently available sug- 
gest an energy dependence in the v e survival probabil- 
ity which is likely inconsistent with the second three- 
neutrino scheme. Taking these two length scales, and 
the results of global flux measurement fits for the solar 
neutrino oscillation interpretation the three differ- 
ent oscillation length and energy scales require three dis- 
parate mass differences, which cannot be accommodated 
in a three-neutrino framework. The CERN e + e _ collider 
LEP measurement of the Z° width indicates the number 
of active neutrinos with masses < mz/2 is 3.00±0.06 0, 
so the results are prima facie evidence for a light sterile 
neutrino species. 

A number of neutrino mass models can provide the 
masses and mixings needed to accommodate all of the 
neutrino oscillation data lEl Ell 51 |23l I2I, ET 
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example, in some string theories, higher-dimensional op- 
erators, suppressed by the powers of the ratio of some 
intermediate mass scale and the string scale, can give 
the light and comparable Dirac and Majorana masses 
necessary for appreciable active-sterile neutrino mixing 
p0[ . In theories with light composite fermions, several 
of the fermions may mix with standard model neutri- 
nos, giving light active and sterile neutrinos UM. In the 
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minimal supersymmetric standard model (MSSM) with 
explicit i?-parity violation, a neutralino can provide the 
required mixing for the atmospheric and solar neutrino 
problems p^ , P3f . This model can also supply a sterile 
or "weaker-than-weakly" interacting particle that has a 
small mixing with one or more active neutrino flavors. 
A low-energy extension of the Standard Model with an 
50(3) gauge group acting as a "shadow sector" may re- 
sult in a neutral heavy lepton [^4| . A model with several 
sterile neutrino dark matter "particles" arises in brane- 
world scenarios. It invokes bulk singlet fermions coupling 
with active neutrinos on our brane j25|. However, many 
models of bulk neutrinos as sterile neutrino dark mat- 
ter must be rather finely tuned and must avoid several 
cosmological constraints 

Although these mass models have been proposed to ac- 
count for the existing neutrino anomalies, many of them 
already contain or can be easily be extended to contain 
additional singlet states which also mix with active neu- 
trinos. As long as the new particles are sufficiently mas- 
sive and have sufficiently small mixings with active neu- 
trinos, they evade terrestrial constraints, and it is in- 
teresting and useful to speculate on their cosmological 
and astrophysical consequences. As we show later in this 
paper, even apparently negligible active-sterile neutrino 
mixing is sufficient to induce the production of sterile 
neutrino dark matter in the early universe. 

A specific model's viability in producing a sterile neu- 
trino dark matter candidate depends on whether the 
mass and mixing properties of the candidate (s) with the 
active neutrinos lie within the range that produces an ap- 
propriate amount of dark matter, and whether the can- 
didate is stable over the lifetime of the universe, does 
not engender conflicts with observationally-inferred pri- 
mordial light element abundances, does not violate CMB 
bounds, and does not contribute excessively to the DE- 
BRA in photons [^7j . Potential constraints on these sce- 
narios also may arise from deleterious effects associated 
with the neutrino physics of core-collapse (Type Ib/c, II) 
supernovae (see, e.g., Refs. f^], [I(J and Sec. |x|). Con- 
straints on massive sterile neutrinos (lOMeV < m s < 
100 MeV) from the SN 1987A signal and BBN were con- 
sidered previously in Ref. |j3lf . 

The existence of massive sterile singlet neutrinos that 
mix with v e has been probed in precision measurements 
of the energy spectrum of positrons in the pion decay 
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., |p2fl , which constrain the mixing matrix element 
2 < 10~ 7 for sterile neutrino masses 50McV < to, < 



130 MeV. These limits may be significantly improved 
in future precision experiments |pq] . Less stringent con- 
straints from peak and kink searches exist for smaller m s . 
In addition, searches for decay of massive v s have yielded 



constraints for \U U 
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|, as well as \U e . 



III. STRUCTURE FORMATION 

Conflicts may have appeared between standard cold 
dark matter theory and simulations and observations of 
large and small scale cosmological structure. Simulations 
predict about 500 small halos (i.e., dwarf galaxies) with 
mass greater than 10 s Mq around a galaxy like the Milky 
Way, but only II candidates are observed near the Milky 
Way [35 1, and only 30 in the local group out to ~1.5 
Mpc [ 36 1 . In other words, simulations of the standard 
A-CDM model may predict more dwarf galaxies than are 
seen (but see Ref. |37|]). Also, even if the dwarf halos 
are dark, their overabundance may hinder galactic disk 
formation [ p8| . 

Another potential problem with CDM simulations is 
the appearance of singularities or "cusps" of high den- 
sity in the cores of halos. The observations of the in- 
nermost profiles of galaxy clusters are ambiguous p9| , 
but rotation curves of the central regions of dark matter- 
dominated galaxies consistently imply low inner densi- 
ties EJJ. Recent TV-body calculations of the nonlinear 
clustering of WDM models have found that enhanced 
collisionless damping can lower halo concentrations, in- 
crease core radii, and produce far fewer low mass satel- 
lites Gil. Also, the observed phase space density in dwarf 
spheroidal galaxies may suggest a primordial velocity dis- 
persion like that of WDM Q. 

A promising new constraint on the nature of dark mat- 
ter may come from study of the Lya forest in the spec- 
trum of high-redshift quasars EJ^, |4| . The structure of 
the Lya forest at high redshifts has been used to con- 
strain the contribution of HDM and, therefore, the mass 
of the active neutrinos |4j| . The constraints presented in 
Ref. Ea] and the Gerstein-Zeldovich-Cowsik-McLelland 
bound [}46| only pertain to fully populated active neu- 
trino seas with a thermal energy spectrum. For ex- 
ample, if an active neutrino has a mass m„ a ~ 1 keV 
(a = /i, t) and the reheating temperature of inflation is 
low (Trh ~ lMeV), then the corresponding active neu- 
trino sea will not be fully populated 47, Q and can be 
WDM EJJ . However, observations of the power spectrum 
of the Lya forest cannot only constrain HDM scenarios, 
but also may be able to test the viability of WDM scenar- 
ios. Since the absolute normalization of the power spec- 
trum is uncertain, the relative presence of power between 
scales near ~ 1 Mpc and near ~ 100 kpc may constrain the 
WDM scenarios that have a relatively large contrast in 
power between these scales |5(|. In addition, the large 
number of high-redshift quasars found by the Sloan Dig- 
ital Sky Survey (SDSS) plj can add considerably to the 
knowledge of the Lya power spectrum, particularly at 
the largest scales. The combination of the galactic power 
spectrum from SDSS and the Lya power spectrum will 
produce even stronger constraints on the behavior of dark 
matter on small scales. 
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IV. PRIMORDIAL LEPTON ASYMMETRY 

The lepton number or asymmetry of a neutrino flavor 
a is defined to be 



(4.1) 



where n Va is the proper number density of neutrino 
species v a , and n 7 = 2£(3)T 3 /7r 2 sa 0.243T 3 is the proper 
number density of photons at temperature T. The lepton 
number(s) of the universe is (are) not well constrained by 
observation. The best limits come from the energy den- 
sity present during BBN and the epoch of decoupling of 
the CMB H, jjH g, H). In fact, the best current bounds 
on the lepton numbers come from the observational limits 
on the 4 He abundance, radiation density present at the 
CMB decoupling, and structure formation considerations 

M M B3: 



-4.1 x 10"^ < L v . < 0.79, 



< 6.0. 



(4.2) 



(4.3) 



The bound on positive L„ c is weaker than that for neg- 
ative L„ e since it is possible to combine the effects of 
neutron-to-proton ratio (n/p) reduction of positive L Vi _ 
with a large L v or L Ut1 which increases the expan- 
sion rate and thus the n/p ratio entering BBN. This 
cancellation could provide a neutrino "degenerate" BBN 
that could replicate not only the primordial 4 He, but 
also the D/H and 7 Li abundances predicted by stan- 
dard BBN. See Ref s. |p2 j, 54fl for a further discussion. As 
stated, the limits (^.2|), (4T3[) depend on an assumption 
of a roughly Fermi-Dirac, low-chemical-potential energy 
spectrum for each neutrino species. They do not ad- 
dress neutrino mass and do not take account of potential 
bounds stemming from closure (i.e., age of the universe) 
considerations. The above limits do not change signifi- 
cantly when recent estimates of observationally-inferred 
primordial abundances are employed In any case, 

it is clear that neutrino number density asymmetries less 
than about 10% of the photon number are easily allowed. 
However, we will see that lepton asymmetries at this 
level, or even several orders of magnitude smaller, could 
have a significant and constraina ble effect if there are 
massive sterile neutrinos (see Sec. VII). 

An issue which arises whenever the lepton number (s) 
L Va differ from the baryon number B (or n = rtf,/n 7 ss 
2.79 x 10~ 8 r2&/i 2 ~ 10~ 10 , where rib is the proper baryon 
number density and fij, is the baryon rest mass clo- 
sure fraction and h is the Hubble parameter in units 
of lOOkms -1 Mpc -1 ) concerns the process of baryoge- 
nesis. For example, electroweak baryogenesis predicts 
the equality of these numbers B — L = 0. Note that 
there also exist simple processes that violate B — L and 
create lepton and/or baryon number either through the 
Affleck-Dine mechanism |56| ] or through non-equilibrium 
decays of a heavy Majorana particle p7f. In Ref. [p8|, 



several natural scenarios of producing a large lepton num- 
ber and the observed small baryon number were investi- 
gated. Furthermore, lepton number could arise sponta- 
neously through matter-enhanced active-sterile neutrino 
transformation at energy scales below that of the baryo- 
genesis epoch |5{| p0|| . Finally, it should be recognized 
that through cancellation of lepton numbers we could 
have B — L = (where L = L Ve + L v ( + L Vt ) , while still 
havin g significant lepton-driven weak potentials [see Eq. 
Q5.12 )] in the early universe. 



V. MATTER-AFFECTED NEUTRINO 
TRANSFORMATION 

A. General framework 

Neutrino mixing phenomena arise from the non- 
coincidence of energy-propagation eigenstate and the 
weak (interaction) eigenstate bases. Eigenstates of neu- 
trino interaction (flavor) include the active neutrinos (v e , 
Vft' v t) w hich are created and destroyed in the standard 
model weak interactions, as well as sterile neutrinos {e.g., 
v s , v' s , v' s ', ...) which do not participate in weak inter- 
actions. Eigenstates of neutrino propagation are states 
of definite mass and energy (or momentum) and evolve 
independently of each other between weak interaction 
vertices. If the bases spanned by these sets of eigen- 
states happen to coincide, then active and sterile neu- 
trinos propagate independently between interactions. In 
general, however, the bases need not coincide, since the 
symmetries of the standard model and its many proposed 
extensions do not require the unitary transformation be- 
tween the bases to be the identity transformation. As 
a result, neutrinos can oscillate, or transform in flavor, 
between interactions. 

This physics applies for any number of neutrino fla- 
vors, including the four (three active plus one sterile) 
which can accommodate the neutrino experiments. In 
the environments we consider, active-active mixing is 
suppressed due to the very similar matter effects for all 
active species. In general, the 4-neutrino evolution may 
have multiple active species mixing with the sterile neu- 
trino, or in multiple-sterile scenarios, there may be mix- 
ing between the sterile neutrinos. In our analysis here, 
we consider the simplifying limit of two neutrino (active- 
sterile) mixing to explore the basic physics of sterile neu- 
trino dark matter production. 

In the case of two-neutrino mixing, the unitary trans- 
formation between the bases can be written as 

\v a ) = cos#|^i) + sin 1 1/2) 
\v s ) = — $m6\v\) + cos0|zaj) 

where \u a ) and \v s ) are active (a = e,/x,T) and ster- 
ile neutrino flavor eigenstates, respectively, and \i>i) and 
^2) are neutrino mass (energy) eigenstates with mass 
eigenvalues m\ and m%, respectively. The vacuum mix- 
ing angle 9 parametrizes the magnitude of the mixing 
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(and, as we shall see, the effective coupling of the sterile 
neutrino in vacuum). We choose all of the neutrino fla- 
vor and mass eigenstates to be eigenstates of momentum 
with eigenvalue p. Then a mass-energy eigenstate \ui) 
{i = 1, 2) develops in time and space with the phase 

(5.1) 



where Ei — y/ m 2 + p 2 is the energy of the eigenstate; 
p = |p| is the magnitude of the proper momentum of the 
species. If the neutrino mass eigenstates are relativistic 
so that Ei 3> rrii, we have B, ~ p + m 2 /2p and x w t, 
yielding the last approximation in Eq. (|5.l|). (In this 



last approximation we suppress the part of the evolution 
operator proportional to the trace, as this simply gives 
an overall common phase to the states.) In our study 
of sterile neutrino production, the sterile neutrinos are 
always relativistic during the epochs in which they were 
produced. 

The difference of the squares of the vacuum neutrino 
mass eigenvalues is, for example, dm 2 = m 2 , — m\. We 
can follow the evolution of a coherently propagating neu- 
trino state Y&v) in either the mass-energy or flavor ba- 
sis. In the flavor basis, a Schrodinger-like equation de- 
scribes how the flavor amplitudes, a a (x) = {v a \$> v (x)) 
with a — e,/i,r and a s (x) = (i/ s \^ ll (x)}, develop with 
time-space coordinate x: 



dx 



ip 



V(x,p) 



V(x,p) - A(p)cos26» A(p)sin20 

A(p)sm28 A(p)cos20~V(x,p) 



r 



(5.2) 



where the first term is proportional to the identity and 
A(p) = 5m 2 /2p. In the context of the early universe it is 
most convenient to take x — t, the Fricdman-Lcmaitrc- 
Robertson- Walker coordinate time (age of the universe) , 
while in supernovae we take x to be position. The weak 
potential V(x,p) represents the effects of neutrino neu- 
tral current and charged current forward scattering on 
particles in the plasma carrying weak charge. In the 
early universe we have V(T,p), but in supernovae, V 
also depends on position x. We suppress the x and T 
dependence of V in the rest of the paper. For a review 
of these issues and neutrino astrophysics, see, e.g., Refs. 

Neutrino mixing can be modified by the presence of 
a finite temperature background and any asymmetry in 
lepton number. The oscillation length is 



I, 



{ A 2 (p) sin 2 26> + [ A (p) cos 29 - V D - V T (p)] 2 } 



-1/2 



(5.3) 



The effective matter-mixing angle is 

A 2 (p) sin 2 26 



sin 29 m = 



A 2 {p) sin 2 26» + [A{p) cos 29 — V D — V T (p)f ' 



(5.4) 



Matter effects have been separated into finite density and 
finite temperature potentials, V D and V T (p). 

The finite density potential V D arises from asymme- 
tries in weakly interacting particles (i.e, nonzero lepton 
numbers), not from non-zero total densities alone. In 
general, the finite density potential is 0, Fiif 



!V2G F [2(n Ve - rip J + {n Vlt - n Vll ) + (n„ T 
V2G F [{n Ve -n Pe ) + 2(n 1/ „ - n p J + {n^ 
V2G F [(n„ e - n Be ) + (n^ - n^) + 2(n Vr 



+ (n e - - n e +) - n n /2] for v e ■ 
rip T )-n n /2] for ^ 

np T )-n n /2] for v T 



(5.5) 



The thermal potential V T arises from finite tempera- of thermally created particles 
ture effects and neutrino forward scattering on the seas _ 

%V2G F Pu 



V T ( P ) = - {{E Va )n Va + {E Da )n Da ) 



'z 
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where n a (n a ) is the proper number density of leptons 
(anti-leptons) of flavor a and (E a ) ((E & )) is the aver- 
age energy of the lepton (antilepton), and n Va (n,p a ) and 
(E Ua ) ({Ep a )) are the proper number density and aver- 
age energy of the neutri nos (antineutrinos) of flavor a. 
The second term in Eq. (5J3) must be included whenever 
the lepton of the same flavor as the active neutrino in 
question is populated. 

In this paper, we will assume that the initial neutrino 
distribution functions are close to Fermi-Dirac black bod- 
ies, which for occupation of differential interval dp have 
the form 



dn„ 



where E(p) — (p 2 



p 2 dp 



T 3 F 2 J e £ (p)/ T -^ + 
/ n 7 \ p 2 dp 



(5.7) 



re 2 ) 1 / 2 , and E(p) « p in the rela- 
tivistic kinematics limit. In this expression i] Va — [i Uoc /T 
is the degeneracy parameter (chemical potential divided 
by temperature) for neutrino species v a and Fzlyj) = 
/ °° x 2 dx/(e x ~ v + 1) is the relativistic Fermi integral of 
order 2 [F 2 (0) = (3/2)C(3)]. The distribution function 
for a neutrino species a is 



f a (p,t) = l/(e E ^ T -^ +1). 



(5.8) 



The last approximation in Eq. (5.7) follows if we take 
the neutrino degeneracy parameter to be zero. This is 
frequently a good approximation for almost all of the 
range of lepton numbers which are interesting for our 
purposes. However, it may not be valid over the bro ader 
range of allowed lepton numbers given in Eqs. (4.2) and 
(4.3). For small lepton numbers r\ Va » 1.46L„ a ; in fact, 
whenever v a has relativistic kinematics, we can write 



4((3) 



(5.9) 



This equation can easily be inverted to find rj Ua {L Va ). 

In the early universe, for temperature ranges where 
the number of degrees of freedom is constant, the time- 
temperature relation is a simple power law and the quan- 
tity e = p/T is a comoving invariant. The differential 
number density for small lepton number in this case is 
then 



dn v 



( n ~i 



e 2 de 



\4C (3) y e e + l 



(5.10) 



B. Quantitative formulation for the early universe 

If one makes the assumption that the only net lepton 
number in the universe is that required for electric charge 



neutrality (i.e., half of the baryon number when there 
are equal numbers of protons and neutrons), then the fi- 
nite density potential usually remains negligible. In some 
cases where Sm 2 < 0, an initially small asymmetry like 
this can be amplified by matter-enhanced active-sterile 
conversion (5^, |6^] . 

The finite density potential in the early universe could 
be dominated by asymmetries in the lepton number, and 
so is often referred to as the "lepton potential." It takes 
the form 



(5.11) 



where we take "+" for a = e and "— " for a = fi, r. Here 
we define the net driving lepton number C a in terms of 
the lepton numbers in each active neutrino species as 



C a = 2L Va 



Lv p 



(5.12) 



with the final sum over the activ e neutrino flavors other 
than v a . Note that in Eq. ( |5.1l| ), the baryon-to-photon 
ratio is 77, not to be confused with neutrino degeneracy 
parameter. 

The total weak potential, V{p, T) = V D {T)+V T (p, T), 
experienced by an active neutrino v a is approximately 



V(p, T) « (40.2 eV) 



f C \f T 



VlO"V VGcV 



Bp 



T 
GeV, 
(5.13) 



As noted above, the thermal potential V T must take into 
account the presence of populated leptons of the same 
flavor. For v e this is required at all temperatures where 
the neutrinos are coupled; for the thermal muon term 
should be included at temperatures T > 20 MeV, where 
the fi is populated; and for v T , the thermal term should 
be included at T > 180 MeV. Therefore, the coefficient 
B takes on the values: 



B 



for T < 20 MeV; 



B 



10.79 eV a = e 
3.02 eV a = (i,T 



(5.14) 



10.79 eV 
3.02 eV 



for 20 MeV < T < 180 MeV; 



10.79 eV a = e,fi,T 



(5.15) 



(5.16) 



for T > 180 MeV. 

The trend of the weak potential experienced by an 
active neutrino species v a is clear. When the quantity 
C ± ?y/4 is sufficiently large and positive (negative for 
Va) the potential will rise with increasing temperature, 
reach a maximum, and then turn over and eventually 



7 



1 



7 ' ' ' 1 1 ' ' 1 


i i i i i i 1 i i i 


i 1 1 1 1 1 




v S3 






Va/ 


\ : 




/ 2 


\ - 




v Si 




r , , , _j_^-rT 





i , ,1 , , j 



T (100 MeV) 

FIG. 1: An example of the temperature evolution of the ac- 
tive and sterile neutrino effective mass-squared, m^. The 
active neutrino is dominated by a positive finite density 
potential at lower temperatures and turns over when the ther- 
mal potential dominates. Resonance occurs at level crossings, 
where the active and sterile m^g tracks intersect. 



become negative at a high temperature where thermal 
terms dominate. The potential as a function of tempera- 
ture for some representative parameters is shown as the 
solid line in Fig. [l| With this behavior it is obvious that 
neutrino mass level crossings in the temperature regime 
not dominated by the thermal terms are possible only 
if the vacuum mass eigenvalues most closely associated 
with some sterile neutrinos are larger than those most 
closely associated with v a . 



VI. NON-EQUILIBRIUM PRODUCTION: THE 
STERILE NEUTRINO BOLTZMANN 
EQUATION 

The Boltzmann equation gives the evolution of the 
phase-space density distribution function f(p, t) for a 
particle species. For sterile neutrinos in the early uni- 
verse, it can be written as |6(| 67 



r t (p' a ,p)f s ( P ,t) [i - f a (p' a ,t)}d 3 P ' a , (6.i) 



where the rate of scattering production of v s correspond- 
ing to a particular channel i [see Eq. (6.8)] is T;. We 
make two approximations regarding the scattering ker- 
nels Ti(p',p): (1) they are isotropic; and (2) they are 
conservative. 



The first of these approximations is completely justi- 
fied given a homogeneous and isotropic universe. The 
second — that an active neutrino scatters into a sterile 
state of the same energy — is made purely to ease the 
computational complexity of following the evolving sys- 
tem of neutrinos. We note, however, that Ref. |38] es- 
timates the effects of relaxing this approximation in the 
context of a semi-analytic calculation and find that it 
makes little qualitative or quantitative difference in the 
results. 

Certainly, however, a better treatment could be ex- 
tended to the scattering kernel. This could, for example, 
change the conditions required for coherence. For exam- 
ple, we note that the effects described in Ref. [|9| may be 
important, especially regarding wher e coherence breaks 
down. We show, however, in Sec. VT] that coherent pro- 



duction of sterile neutrinos at MSW resonances is not 
important in our favored parameter regime. 

In our discussion, we take fi (p, t) to be the momentum- 
and time-dependent distribution function of active (a = 
e, fi, t) and sterile (s) neutrinos. The distribution func- 
tions f a are given by Eq. (Ol) . The second left hand side 
term arises from the redshifting of the distribution. 

However, Eq. ( |3.1| ) is semi-classical in its evolution of 
the neutrino distributions. To exactly follow the full 
quantum development, including such quantum effects as 
damping and off-diagonal contributions to the neutrino 
Hamiltonian in matter |70|| , one must follow the time 
evolution of the density matrix |7l|, |72|, |7^, Q . One can 
approximate the effects of quantum damping through a 
damped conversion rate T(v a — > v s ;p,t). Fermi blocking 
effects are taken to be negligible, leaving the Boltzmann 
equation as 



d_ 

df 



Hp—f s {p,t) 



T{v a -> u s ;p,t) [f a (p,t) - f.fat)] . (6.2) 



The reverse term is important for the L ^ case, where 
f a and f s can become comparable for certain regions of 
the neutrino momentum distributions. 

The conversion rate to sterile neutrinos is just the 
product of half of the total interaction rate, T a , of the 
neutrinos with the plasma and the probability that an 
active neutrino has transformed to a sterile: 

T(v a -> V s \p, t) W -^-{PmiVa -> V s ]P,t)), (6.3) 

This probability P m depends on the amplitude of the 
matter mixing angle and the quantum damping rate 
D(p) — Tg(p)/2 \D(p) = T a (p)/2] for neutrinos [antincu- 
trinosl ifM vQ, Wft, U^l, 
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(P m (va -> v s ;p, t)) 
(P m (v a -> v s ;p, t)) 
The full Boltzmann equation then is 



1 A(p) 2 sin z 26 

2 A(p) 2 sin 2 26 + D 2 (p) + [A(p) cos 26* — V D — V T (p)} 2 

1 A(p) 2 sin 2 26» 

2 A(p) 2 sin 2 26 + D 2 {p) + [A(p) cos 26> + v" 5 - V T (p)} 2 



(P m (v a -> ^s,£in + T"))r [fa{p,t) ~ f s (j>, t)] 



' 26,, 



T a (j>)lr. 



T a {p)A 2 (p) sin 2 26> 



(6.4) 
(6.5) 



(6.6) 



4 A 2 (p) sin 2 26> + D 2 {p) + [A(p) cos26» - F L - ^ T (p)]' 



[/„(?,*)- /.(p,t)] 



r 



where {1 + [r a (p)Z m /2] 2 } -1 is the damping factor. T here 
are analogous equations for Eqs. (6.1)-(6.3) and Eq. ( |6.6| ) 
for antineutrinos. 

Previous calculations have approached the solution of 
this equation for the case of negligible lepton number 
(L w 0) analytically and have been restricted to the short 
epoch just prior to BBN where the Hubble expansion rate 
H (evolution of the scale factor) and time-temperature 
relations arc simple power laws (the temperature is pro- 
portional to the inverse scale factor) ^38|. Such ap- 
proximations allow the reduction of the left-hand side of 
Eq. ( |6.6| ) to a single term. With the approximations that 
lepton number C is always negligible, that the thermal 
term V T is not modified by population of leptons, that 
the interaction rate is not enhanced due to population 
of scatterers, that quantum damping is never important, 
and that the reverse rates (y s — > v a ) are always negli- 
gible, then the right hand side can also be considerably 
simplified, and the solution for the sterile neutrino dark 
matter abundance is reduced to a simple integral. It is 
obvious, however, that many if not all of these approxi- 
mations are eventually invalid over at least some of the 
parameter range of interest for sterile neutrino dark mat- 
ter. 

For example, in order to probe cases where sterile neu- 
trino dark matter production may lie above the QCD 
transition (T > 100 MeV), and in order to improve the 
accuracy of the predicted dark matter contribution, we 
must extend our calculation to epochs where the scale 
factor-temperature and time-temperature relations are 
not simple power laws (see the Appendix for details). 
The interactions contributing to T a (p) which produce 
(and remove) sterile neutrinos are 



v a + v [3 
v a + l± 
v a +q 
v n + v n 



i+ + r. 



(6.7) 



Here, the v a represent either neutrinos or anti-neutrinos, 
as appropriate, and q and / are any populated quark and 
charged lepton flavors. The total interaction rate at tem- 
peratures 1 MeV < T < 20 MeV due to interactions of 
the neutrinos among themselves and the pairs is 



r Q (p) 



1.27 G|pT 4 , 
0.92 G 2 pT 4 , 



a 
a 



e, 

fi,r. 



(6.8) 



At higher temperatures, other leptons and quarks are 
populated and contribute to the neutrino interaction 
rate. In our calculations, we have included the enhance- 
ment of the interaction rate due to the presence of these 
new particles in the plasma at high temperatures. In par- 
ticular, a significant increase to the scattering rate results 
at temperatures above the QCD scale. Interestingly, the 
results of the production of sterile neutrino dark matter 
therefore depend on the temperature of the QCD transi- 
tion, where the quark-antiquark pairs annihilate and are 
incorporated into color singlets. 

The high scattering rate characteristic of the environ- 
ment of the early universe not only serves to populate the 
sterile neutrino sea, but can also suppress the production 
of sterile neutrinos when the matter oscillation length 
is large compared to the mean free path of the neutri- 
nos. When the oscillation length is much larger than the 
mean free path, the probability that an active neutrino 
has transformed into a sterile state becomes very small. 
It can be shown that such scattering can force a quan- 
tum system to not evolve from the initial state |75|, [79] . 
Essentially, each scattering event resets the phase of the 
developing neutrino state l^^)- This is the so-called 
quantum Zeno effect. 

Consider the right hand side collision term in Eq. (|6.2|) . 
The sterile neutrinos are initially not in thermal equi- 
librium, and therefore the reverse processes are initially 
unimportant. They will, however, become more impor- 
tant as the sterile neutrino sea is populated. Because 
the sterile neutrinos are never in equilibrium, the usual 
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simplifying principle of steady-state equilibrium fl67| , pQ ] 
cannot be made. Therefore, to calculate the production 
of sterile species, we must start with t he B oltzmann equa- 
tion in its "unintegrated" form Eq. (3.6). 

We have directly solved the Boltzmann equation nu- 
merically for the L w case {e.g., L w 1CP 10 ), and found 
the amount of sterile neutrino dark matter produced for 
a broad range of sterile neutrino mass and mixing angles 
with each of the three active neutrino flavors. In treating 
the redshift of the sterile distribution, we have greatly ac- 
celerated the numerical calculation by following the red- 
shift of the sterile neutrino distribution function f s by 
redshifting momenta as p oc i -1 / 2 , which is always true 
in a radiation dominated universe (a necessary condition 
for BBN and the CMB). This numerical method allows 
us to correctly follow the redshifting of the sterile neu- 
trino without the second ter m o f the Liouville operator 
in the Boltzmann equation ( |6.6[ ), while active neutrinos 
may be reheated via annihilation of other species. 

The contribution to the closure fraction of the universe, 
meeting all constraints (see Sees. IX and |x|) is shown in 
Fig. U for v e ^ v s and Fig. || for v T ^ v s . A general fit 
to our results for nonresonant production is 



n„ h 2 



0.3 



/ sin 2 


26» \ 


/ m s \ 


1 io- 


10 J 


UOOkeV/ 



(6.9) 
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FIG. 2: Regions of Q, Vs h 2 produced by resonant and nonres- 
onant z^ e <-» v s neutrino conversions for selected net lepton 
number L, after applying all constraints (see Sees. |k| and 
[xj). Regions of parameter space disfavored by supernova core 
collapse considerations are shown with vertical stripes. 



The maximum rate of sterile neutrino production occurs 
at temperature pOl plf 



133 MeV 



(l 



1/3 



keV/ 



(6.10) 



The closure fraction contribution for ^ v s is nearly 
identical to v T # v s . The particular flavor (s) of the ac- 
tive neutrino with which the sterile neutrino mixes does 
determine to some extent the ultimate closure fraction. 
However, this flavor dep enden ce is negligible for larger v s 
masses since, from Eq. ( |3.1C ), for sterile neutrinos with 
m s > 2keV, T max > 180 MeV. At these temperatures, 
e^, fj, and even (due to the high entropy of the 
universe) are populated significantly, so that the ther- 
mal potentials are identical for all flavors. Therefore, the 
mass fractions produced for all flavors with sterile neu- 
trino masses m s > 2keV are very similar. 

The similarity of the results of production for v e , z^, 
and v T mixing with a sterile neutrino is at odds with 
the calculation of Ref . |68| . There are several differences 
between our treatment and that of Ref. J6^] that can 
account for the disparity in results. First, the primary 
cause for discrepancy for the v T ^ v a neutrino mixing 
case is likely to lie in the fact that r leptons are signifi- 
cantly populated at temperatures one-tenth of their mass 
(T > 177 MeV) due to the high-entropy of th e universe. 
This modifies the thermal term in Eq. (5.6) in an im- 
portant way, not included in Ref. |38|. In our work, we 
follow the number and energy density of tau leptons ex- 
plicitly in our numerical evolution. One can see from Eq. 
(3.10) that for masses of sterile neutrinos greater than 



about 2 keV, production occurs at a temperature where 
the t lepton is significantly populated. Second, the pop- 
ulation of massive species of scatterers (p and r leptons 
and u, d, s, c, b quarks) which enhance the scattering rate 
is not included in Ref. [ |68[ , but are included in our cal- 
culations. Third, effects of re- heating on the dilution 
of sterile neutrino dark matter are treated only approxi- 
mately in Ref. [Q. We explicitly include reheating and 
dilution in our calculation through our treatment of the 
time-temperature evolution of active and sterile neutri- 
nos, as described in the Appendix. 



VII. INITIAL LEPTON NUMBER AND 
RESONANT PRODUCTION 

In this section, we examine the resonant production of 
sterile neutrino dark matter for a range of initial lepton 
numbers. We will extend the treatment of this issue given 
in Shi and Fuller |ll| considering the effects of particle 
annihilation and reheating and enhanced scattering rates 
on coherent and incoherent evolution through MSW reso- 
nances. This allows us to consider resonant v a ^ v s con- 
version at high temperature epochs in the early universe, 
up to the electroweak transition regime at T 100 GeV. 

Resonance, or mass level crossings as seen in Fig. [l], arc 
characterized by maximal effective matter mixing angles 
(0m)res = 7r/4, where the v a — > v s conversion rate is, 
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consequently, enhanced. The resonance condition is 

A(p) cos 29-V L - V T (p) = 

'£±77/41 / T 



A(p) cos 26»- (40.2 eV) 



10" 



GcV 



VGeV 

2 



= 



(tSv)'» s2 ^ 8 - 03 «{^}(too^v)' 



v_zi_y 

VkcV/ VlOOMeV/ 



(7.1) 



Choose one of the horizontal dotted lines laying be- 
low the peak of the solid line in Fig. [j]. This indicates 
the "mass track" for a sterile neutrino species. This 
is, of course, simply the vacuum mass-squared value 
(m 2 s ) s — m 2 , and is flat and independent of tempera- 
ture. The effective mass-squared track for an active neu- 
trino v a could be as shown in Fig. |l|. Mass level crossings 
(resonances) in the v a ^ v s system can occur when m 2 
lies below the peak value of (m 2 s ) Ua . This peak will oc- 



cur at temperature 



Tpeak 



1/2 



1/2 



/4^2C(3)\ 
G F (GeV) 5 



-1/2 



{£±7?/4} 



1/2 



B 



(2.98GeV)e- 



1/2 



-1/2 



C 



10" 



1/2 



(7.2) 



where To is the high temperature at which the effective 
mass-squared track crosses zero, and where in the last nu- 
merical expression we have assumed that T > 180 MeV. 

Note that the peak value of effective mass-squared, or 
equivalently, the largest value of Sm 2 cos 28 for which a 
level crossing can occur is, 



m 2 cS (Tpeak) = (£m 2 cos20) 



cff(PEAK) 



/ WgqgA {£±r;/4} 3 

I 3vr 2 J e 



(7.3) 

G f (GeV) 10 



D 2 



For example, we can show that the largest u s mass which 
can have a resonance with a v e is roughly 



im s ) 



406 keV/ L 



PEAK 



.1/2 



Go 



3/2 



(7.4) 



As the universe expands, an active neutrino species 
v a will encounter two resonances with a sterile neutrino 
species, as long as this sterile s peci es has m 2 less than 
the peak value of m 2 s in Eq. (7.3). At some epochs, 



both resonances may be present in a given active neu- 
trino spectrum since the resonance condition, Eq. ( |7.l| ), 
has more than one zero. This behavior is shown for a 
particular case in Fig. ||. 

What is the effect of resonance on sterile neutrino pro- 
duction in the early universe? The answer to this ques- 
tion depends on the rate of incoherent production at the 
resonance, whether the neutrinos are coherent as they 
pass through resonance and, if they are coherent, the de- 
gree of adiabaticity characterizing the evolution of the 
neutrino flavor amplitudes through the resonance. There 
are two processes by which resonance affects neutrino 
conversion in the early universe: (1) enhanced incoherent 
conversion of those neutrinos at the resonance and (2) co- 
herent MSW transformation of neutrinos passing through 
the resonance. The rate of incoherent production can be 
calculated numerically through the semi-classical Boltz- 
mann evolution, Eq. (16^). 

The statistical formulation of both coherent and in- 
coherent production can be described by the time evo- 
lution of the density matrix for the two neutrino states 
@ [Fl], [72| [73|, |74|, fF§. We show below that the masses 
and mixing angles of interest for sterile neutrino dark 
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FIG. 4: The effective matter mixing sin 28 m is shown for a 
case with two resonances (solid line). For this case, m a = 
IkeV, sin 2 26 = 1CT 9 , L = 6 x 1(T 4 , and T = 130 MeV. The 
dashed line is the active neutrino energy distribution. 



matter considered here give coherence across the rele- 
vant resonance width, but usually imply that the neu- 
trino amplitude evolution through resonance is nonadia- 
batic. (Adiabatic evolution at lower mixing angle could 
occur at epochs where the entropy is being transferred 
from annihilating particles and where, as a consequence, 
the temperature and density do not change rapidly or are 
constant with time, e.g., the QCD transition fr&| .) 

Though the expansion rate of the universe scales as 
T 2 , the active neutrino scattering rate scales as G|T 5 . 
For neutrino flavor evolution through resonances, the co- 
herence condition will be met only when the inverse of 
the scattering rate (the active neutrino mean free path) 
is much larger than the resonance width or 



T < 



4tt 3 \ 



1/6 



9 



1/6 



140 MeV- 



(,g/100) 1/6 



sin 2 2#/nr 10 ) 1/6 ' 



(7.5) 



(Here, m p i « 1.22 x 10 22 MeV is the Planck mass.) 
Therefore, for the mixing angles relevant here (sin 2 28 < 
1CP 10 ), at temperatures below 140 MeV, the resonance 
could in principle drive coherent v a v s transformation. 
The efficiency of MSW conversion relies on the adiabatic- 
ity of the evolution through the resonance region. 

The width of the resonance is the product of the 
local density scale height of weak charges and tan 29. 
In turn, the biggest share of the density scale height 
is determined by the expansion rate of the universe, 
H w (87r 3 /90) 1/2 3 1/2 T 2 /m p i, where the statistical weight 
in relativistic particles, g, has contributions from both 
bosons <7b and fermions gf. 



S = £(ff&)i + 7/8X>/)* 



(7.6) 



In fact, once neutrino flavor transformation begins, the 
inherent nonlinearity of this process can have a sizable, 



even dominant effect on the density scale height. Ignor- 
ing this, the resonance width expressed in time is St w 
(2/3)t(tan20) = (l/tyH- 1 tan 2(9, where t « (l/2)fl"- x 
is the age of the universe at an epoch with temperature 
T in radiation dominated conditions. [The particle hori- 
zon is -ff -1 , so that, absent large scale neutrino flavor 
transformation, the resonance width is a constant frac- 
tion (1/3) tan 29 of the horizon scale.] 

The effective matter mixing angle for the oscillation 
channel v a ^ v s at an epoch with temperature T is 



sin 2 2dr. 



[1 - 2eTV/ (Sm 2 cos 20) 
tan 2 29 



(7.7) 



At resonance, sin 29 m = 1, so that one resonance width 
off resonance this effective mixing will have fallen to 



' 26m 



1/2. Clearly, the change in effective weak 



potential over this interval is SV 



! sin20/(2eT res ), 



where T rcs is the resonance temperature for neutrino 
spectral parameter e. From this it can be seen that 
(SV/V) = tan 29 and it follows that the resonance 
width is 



1 dV 
V ' ~dt 



tan 26. 



(7.8) 



The neutrino oscillation length at resonance is 1™ S = 
(4ireT les ) / (Sm 2 sin 29) = 2ir/5V, and the adiabaticity pa- 
rameter is proportional to the ratio 5t/l™ s , and is defined 



St 



7 = 2tt— *(5Vy 




(7.9) 



When there are many oscillation lengths within the reso- 
nance width (i.e., when 7> 1) neutrino flavor evolution 
will be adiabatic and we can then expect efficient flavor 
conversion at resonance. [Given the small vacuum mix- 
ing angles relevant for this work and consequent small 
widths, the Landau-Zener jump probability |83|[ with the 
form Plz ~ exp (—tt^/2) gives an adequate gauge of the 
probability of v a — > v s conversion at r eson ance, 1 — Plz-] 



The physical interpretation of Eq. (7.9) is straightfor- 
ward. The adiabaticity parameter is proportional to the 
energy width of the resonance divided by the "sweep 
rate," de/dt, of the resonance energy through the neu- 
trino distribution function. The resonance sweep rate 
is determined mostly by the expansion rate of the uni- 
verse (an inverse gravitational time scale), but the rate 
of change of lepton number C as v a — > v s proceeds can 
become important, even paramount as lepton number is 
used up and C — > 0. 

The resonance energy width scaled by temperature is 



Se w 5V\de/dV\ 



e res tan 29, where e res is the reso- 



nant value of the neutrino spectral parameter at T Tes . If 
we confine our discussion to resonances on the low tem- 
perature side of Tpeak, where the thermal terms in the 
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potential can be neglected, then 

Sm 2 cos 29 
£rcs * (4V2C(3)/7r 2 )G F T4£ 



(7.10) 
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As the universe expands and cools with time, and for 
a given 5m 2 , the resonance will sweep through the v a 
energy distribution function from low to high neutrino 
spectral parameter e. In this same limit of resonance 
below Tpeak, the sweep rate is 



dt \ AHC 



(7.11) 



where C is the time rate of change of the lepton num- 
ber resulting from neutrino flavor conversion. Since the 
expansion rate scales as H ~ T 2 , the prospects for adia- 
baticity are better at lower temperatures and later epochs 
in the early universe, all other parameters being the 
same. 

Prom Eqs. Q, ( gig ), and ( f7ll| ), we can estimate 
that at resonance the degree of adiabaticity is 



Sm 2 sin 2 26* ^ T . t 



2eT res cos26> \ T C 
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(7.12) 
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where in the second equality we assume the standard 
radiation-dominated conditions and expansion rate, and 
where in the final equality we have employed the approx- 
imation dm 2 w m 2 , valid when m s 3> Tn Uct , and where 
we have assumed that the vacuum mixing angle is small. 
Here we see that for the mixing angles allowed by our 
constraints, sin 2 29 < 10 _9 (3 x 10~ 10 ) for v^, v T {v e ) mix- 
ing with sterile neutrinos, and masses m s > 1 keV, the 
resonance is not adiabatic. 

We conclude that the main effect of resonance is en- 
hancement of scattering-induced incoherent conversion of 
neutrinos with energies in the resonant region. Therefore, 
t he f ormulation of the semi-classical Boltzmann Equation 
(6.6) is appropriate for calculating the total production 
of sterile neutrinos in the early universe. 

The results of our numerical calculations can be seen 
in Fig. ^ for v e ^ v s and in Fig. [| for v T ^ v s for 
the cases where initially L — 0.001,0.01,0.1. The calcu- 
lation includes both nonresonant scattering production 
and matter-enhanced (resonant) production. Examples 




g = p/T 

FIG. 5: The sterile neutrino distribution for four cases of res- 
onant and non-resonant v e <-> v s , as described in the text. 
The dotted line is a normalized active neutrino spectrum. 
The thick-solid, dashed, dot-dashed, and thin-solid lines cor- 
respond to cases (l)-(4), respectively. The inset shows a mag- 
nified view of the low momenta range of the distributions. 



of the resulting sterile neutrino energy spectra are shown 
in Fig. |5|. Resonantly produced sterile neutrinos tend to 
have energy spectra appreciably populated only at the 
low e end. This results from the resonant energy start- 
ing at the lowest momenta and m oving through higher 
momenta neutrinos [see Eq. (7.10)] as the universe cools 
and lepton number is depleted through conversion into a 
sterile neutrino population. 

Figure |^ shows the resulting spectrum for four sample 
cases of sterile neutrino dark matter production: 



(1) m s = 0.8 keV, sin 2 29 = 10~ 6 , i init = 0.01, resulting 
in n Us h 2 = 0.25 and (p/T) = 2.9; 

(2) m s = 1 keV, sin 2 29 = 10~ 7 , L init = 0.01, resulting 
in n u h 2 = 0.13 and (p/T) = 1.8; 

(3) m s = 1 keV, sin 2 29 = 10~ 8 , L init = 0.01, resulting 
in Vl u h 2 = 0.10 and (p/T) = 2.0; 

(4) m s = 10 keV, sin 2 29 = 10~ 8 , L inlt = 0.001, resulting 
in n u h 2 = 0.57 and (p/T) = 2.3. 

A particularly interesting case is (4), where the resonance 
passes through the distribution during the QCD transi- 
tion, where the disappearance of degrees of freedom heats 
the photon and neutrino plasma, forcing the universe to 
cool more slowly. For this period, the resonance moves 
much more slowly through the spectrum and is conse- 
quently more efficient in v a — > v s conversion through 
that region of the neutrino energy spectra. This pro- 
duces a "spike" in the sterile neutrino distribution (see 
Ref. M). 
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VIII. COLLISIONLESS DAMPING SCALE: 
HOT, WARM OR COLD NEUTRINOS 

Given the sterile neutrino energy spectrum and num- 
ber density, the transfer function for dark matter models 
can be calculated. Here, we instead give a rough guide 
based on the free streaming length at matter-radiation 
equality Xfs- Structures smaller than Xfs & r e damped. 
This is because at early epochs where sterile or active 
neutrino species are relativistic, they can freely flow out 
of the regions of sizes smaller than Xps (very roughly 
the horizon size at the epoch where the neutrinos revert 
to nonrelativistic kinematics). Previous numerical work 
has shown that the free streaming scale is approximately 

The mass contained within the free streaming length is 
then 

M ra ,, 6xl O"M (^)(iM) 3 ((|§) 3 , 

(8.2) 
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fi„ h 2 : CD 0.1-0.5 □ o.oi-o.i I I 0.001-0.01 



FIG. 6: The mass within the sterile neutrino free streaming 
length at matter-radiation equality (Mfs) in the v T <-» v s , 
L = 0.01 case. Regions of parameter space disfavored by 
supernova core collapse considerations are shown with vertical 
stripes. 



where fi m is the contribution of all "matter" to closure. 
These values can give a guide as to the collisionless damp- 
ing scale of sterile neutrino dark matter. In Fig. we 
show contours of Mps for the v T ^ v s sterile neutrino 
production channel with L = 0.01. In Fig. ^, we assume 
that the universe is critically closed (Vl m = 1). Therefore, 
the contours are consistent with sterile neutrinos being 
the major constituent of dark matter near the dark gray 
regions where fi„ h 2 = 0.1 — 0.5. 

In Figs. H and ||, we label the parameter regions corre- 
sponding to M FS > 1O 14 M , 10 5 M Q < M FS < 1O 14 M , 
and M FS < 10 5 M Q as HDM, WDM, and CDM, respec- 
tively. These definitions are somewhat arbitrary and here 
serve only as guides. Note that sterile neutrino mass by 
itself does not completely determine the free streaming 
(or collisionless damping) scale, since the inherent v s en- 
ergy spectrum at production also helps to determine the 
v s kinematics at a given epoch. This latter effect is es- 
pecially pronounced for nonthermal energy spectra, and 
is responsible, e.g., for the non-flat collisionless damping 
mass scale lines in Fig. ^. 



A. Diffuse extragalactic background radiation 

The decay rate of a massive sterile neutrino with vac- 
uum mixing angle into lighter active neutrinos is pq] 



.7 x lO^s" 1 
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We have also included in our calculations the contri- 
butions to r„ s from visible and hadronic decays esti- 
mated from the partial decay widths of the Z° boson . 
The rate of the corresponding radiative decay branch is 
smaller by a factor of 27a/87r pq]: 
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IX. COSMOLOGICAL CONSTRAINTS ON 
STERILE NEUTRINOS 



That these sterile neutrinos can decay stems from the 
fact that they are not truly "sterile," but have an effective 
interaction strength ~ sin 2 ^ 29 G| = sin 2 



26» Gf , where 



the last equality is valid in vacuum (i.e., at late epochs 
in the universe). Obviously, if a particle is going to be a 
dark matter candidate, it must have a lifetime r at least 
of order the age of the universe (itoday ^ lOGyr). Re- 
gions of m s vs sin 26 parameter space where r < t to d 
and not otherwise constrained are so labeled in Figs, 
and R These parameters, corresponding to generally 
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high sterile neutrino masses, guarantee that these par- 
ticles decay away into lighter particles long before they 
clump in gravitational potential wells. 

The radiative decay channel for sterile neutrinos can 
provide a constraint or a possible detection mode for 
some regions of the m s — sin 2 20 parameter space. In 
particular, radiative decays of sterile neutrinos occuring 
between CMB decoupling and today could produce an 
appreciable flux of photons with energies comparable to 
m s . In fact, there are firm upper limits on the flux from 
a diffuse photon component. For example, the total dif- 
ferential energy flux per unit solid angle for a DEBRA 
component is p7l p7|, [S7| 



dF/dVL < (lMeV/E) cm" 2 sr-V 



(9.3) 



As technology has progressed this limit on the true dif- 
fuse background has come down as distinct x-ray sources 
are resolved and their fluxes are removed from the count 
p8| , psfl . At present the Chandra X-ray Observatory 
threatens to resolve a considerable fraction of the ob- 
served x-ray background (in the 0.5 — 8keV band) into 
point and extended sources, primarily active galaxies, 
QSO's, and clusters. 

Clearly, sterile neutrinos with radiative decay rates 
greater than an inverse Hubble time (H^ 1 w 3.09 x 
10 17 h~ 1 s « 9.78h~ 1 GyT) will tend to contribute back- 
ground photons before the sterile neutrinos fall into po- 
tential wells and form structure. These photons will then 
produce a DEBRA contribution which must not exceed 
the overall limit in Eq. (9.3). Parameter regions violating 
this bound are labeled DEBRA in Figs. and ||. 

What about sterile neutrinos with much smaller radia- 
tive widths? These steriles could be decaying at more 
recent epochs, even today. However, if these are the 
dark matter (CDM or WDM), then they are not dif- 
fuse, but are strongly clustered. In Figs. and || the 
regions with widely spaced vertical lines correspond to 
sterile neutrino mass and mixing properties that would 
give a DEBRA component in excess of the limit ( |9.3| ), 
if these steriles were distributed diffusely. Since most of 
the decay photons will be produced when the dark mat- 
ter is in structure rather than diffuse, this region does 
not as yet constitute a constraint. Rather, it serves to 
define parameters that could give interesting x-ray fluxes 
in the gravitational potential wells of clusters of galaxies 
or other structures, depending on redshift, the cosmolog- 
ical parameters, and transfer functions and collisionless 
damping scales of the sterile neutrinos. Improved obser- 
vations and models could lead to these regions turning 
into true constraints and may result in more stringent 
constraints or even lead to detection Q . 

As an example, consider a large cluster of galaxies with 
dark matter mass M « 1O 15 M0. The sterile neutrino 
decay luminosity in photons is 



Lk 7 x 10 36 ergs 
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The observed x-ray luminosity for such clusters in the 
1 - 10 keV band is L clustcr - 10 45 erg s" 1 |l). Clearly, for 
sterile neutrinos in our dark matter parameter space with 
masses of 20 — 30 keV and sin 26 w 10~ u , the predicted 
luminosities may be comparable to those observed. The 
energy spectra may be different, however, and this could 
be an avenue for constraint. 



B. Cosmic microwave background 

Another constraint stems from the increase in energy 
density in relativistic particles due to massive sterile 
neutrino decay prior to cosmic microwave background 
(CMB) decouphng |§. The BOOMERanG/MAXIMA 
observations |93j currently limit the effective number of 
neutrinos (N„) at decoupling to iV„(CMB) < 13 at the 
95% confidence level J55|. Measurements to higher mul- 
tipole moments by the MAP and Planck probes will be 
able to further limit the relativistic energy present at 
decoupling. MAP may constrain TV,, (CMB) < 3.9, and 
the Planck mission could reach a limit ^(CMB) < 3.05 
p4j . The increase in energy density due to sterile neu- 
trino decays was found by calculating the energy pro- 
duced by decays between active neutrino decoupling 
(T - 1 MeV) and photon decoupling (T ~ 0.26 eV). The 
BOOMERanG/MAXIMA constraint is shown in Figs. @ 
and U and the potentially stringent future "constraints" 
by MAP and Planck are also shown. 



C. Big bang nucleosynthesis 

The energy density in the sterile neutrino sea at weak 
freeze-out just prior to primordial nucleosynthesis (T sa 
0.7 MeV) must not be too large. The energy density con- 
tribution of a sterile neutrino is often described as the 
fraction of the energy density in a fully populated (ther- 
mal) single neutrino-antineutrino species plus the energy 
density in the active neutrinos, N v . Depending on one's 
particular adoption of observationally inferred primordial 
abundances, one can arrive at limits between N v < 3.2 
and N v < 4 |95|, |9(| ||) . We have calculated the energy 
density contributed by the sterile neutrino at BBN, and 
its limits are shown in Figs. and 



D. Li and D photoproduction 

Another constraint arises from photoproduction of 
deuterium (D) and 6 Li stemming from the decay of mas- 
sive sterile neutrinos after BBN [Q. Energetic cascades 
dissociate 4 He into excessive amounts of D, which is 
bounded observationally p9[ . Also, energetic 3 H and 
3 He produced in the cascades can synthesize 6 Li through 
3 H( 3 He) + 4 He Li + n(p), o verp roducing 6 Li, which 
is also bounded observationally [10C]. 
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These constraints lie in regions of parameter space 
which are almost identical to those bounded by CMB 
considerations. The 6 Li limit can be seen to extend to 
the small region just to the left of the CMB constraints 
in Figs. ^ and || (closely hatched area, hatching oriented 
from the lower right to the upper left). It should be 
noted that the 6 Li/D photoproduction constraints are 
from current observation, while the potential constraints 
from MAP and Planck of the energy density present at 
CMB decoupling are a future possibility that can corrob- 
orate the photoproduction constraints. 
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FIG. 7: Contours of closure fraction are displayed as thick 
lines labelled with Q, Va h 2 for negligible lepton asymmetry 
(L ps 1CP 10 ) for the three neutrino flavors mixing with a mas- 
sive sterile neutrino. Shown are constraints from the energy 
density present at CMB decoupling, with the dark gray region 
corresponding to BOOMERanG/MAXIMA's present limits. 
The medium gray region corresponds to MAP's predicted 
future constraint, and the light gray region is the Planck 
mission's potential constraint region. The constraints from 
DEBRA, and BBN are so labelled. The region of potential 
constraints from supernovae (SN) is also shown as closely 
spaced vertical lines. The region of widely spaced vertical 
lines is where radiative decays today may give detectable X- 
ray signatures. Regions filled with horizontal lines are where 
Q.„ s h 2 > 0.5 and are inconsistent with the observed age of 
the universe. Sterile neutrinos with parameters in the region 
labelled by "r < t t oday" decay without constrainable effects 
and make no contribution to the present matter density. 




FIG. 8: Contours of closure fraction are displayed as thick 
lines labelled with h 2 for L = 0.001, 0.01, 0.1. Constraints 
are shown as in Fig. M. 



X. ACTIVE-STERILE NEUTRINO 
TRANSFORMATION IN CORE-COLLAPSE 
SUPERNOVAE 

In the previous sections, we have calculated the ranges 
of sterile neutrino masses and mixing angles (and pri- 
mordial lepton asymmetries) for which sterile neutrinos 
can account for some or all of the dark matter. In this 
section we describe some of the implications for core- 
collapse (Type Ib/c, II) supernovae of active-sterile neu- 
trino transformation with parameters in these regions. 
Neutrinos play a dominant, pervasive role in core-collapse 
supernovae, so the stakes are high whenever we introduce 
non-standard neutrino physics like neutrino mixing. 

In broad brush the effects are simple: too much neu- 
trino conversion in a supernova results in too much en- 
ergy loss to sterile neutrinos, manifestly in conflict with 
observation (the mixing angles considered here are suffi- 
ciently small that the sterile neutrinos are not trapped in 
the core). Requiring the hemorrhaging in sterile neutri- 
nos not to be too great then places bounds on neutrino 
mixing parameters. 

In fine detail, however, the procedure is not so sim- 
ple, since the coupled supernova neutrino transport- 
transformation problem is highly nonlinear and the dy- 
namics very difficult to treat analytically or numerically. 
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) to describe 



Confounding the matter is our relatively poor under- 
standing of the physics of the core-collapse phenomenon 
itself. 

Nevertheless, we have boldly attemp ted (as have sev- 

" BUT 



eral workers before us 

some of the salient effects of neutrino transformation in a 
supernova and have extracted some rather conservative 
limits on mixing parameters, conservative, at least, for 
arguing the case for sterile neutrino dark matter. Our 
"limits" are by no means final, and much future work 
remains to be done in order ascertain the viability of 
sterile neutrino dark matter. In what follows, we give a 
brief biography of a core-collapse supernova, describe the 
relevant neutrino transport and transformation physics, 
and indicate how we obtained the supernova "constraint" 
regions in Figs, g g g, g, and g. 

Core-collapse supernovae are the death throes of mas- 
sive stars. A star of mass > lOAf© dies when its w 1.4M© 
iron core undergoes gravitational collapse to a neutron 
star, an event releasing « 99% of the « 10 53 erg gravi- 
tational binding energy of the neutron star in all active 
species of neutrinos (see Refs. [101, 102 for an introduc- 
tion). In the standard picture this event lasts ~ 10- 



15s 

and can be divided roughly into three phases: the infall, 
shock reheating, and r-process epochs. 

During the infall epoch, the iron core collapses on a 
near gravitational time scale of ~ 1 s to a dense proto- 
neutron star. As the core density rises, the forward reac- 
tion in 



y 



(10.1) 



neutronizes the core until the neutrino trapping density 
10 12 g cm -3 is reached, a condition which al- 



P 



10 1 



lows the forward and reverse reactions in Eq. (10.1) to 
achieve /^-equilibrium. This equilibrium is established 
with relativistically degenerate electrons (fi e ~ 25 — 
220 MeV) and electron neutr inos (fi Ue w 10 - 170MeV). 

The quotes in Eq. ( 10.1 ) refer to free and bound 
baryons. Most of the baryons are bound in large nuclei, 
since the infall-collapse phase is characterized by tem- 
peratures T rs 1 — 3 MeV and a low entropy-per-baryon 
s w 1.5 (in units where Boltzmann's constant is unity). 

The infall epoch ends when the core density reaches 
and exceeds the saturation density of nuclear matter. 
The inner core bounces, yielding an outward-propagating 
shock wave at its boundary with the outer core. The 
shock loses much of its energy dissociating the nuclei 
in the outer core and mantle of the pre-supernova star 
and eventually stalls some 500 km from the center of the 
core, far short of generating a supernova explosion. The 
post-bounce core is a hot, dense proto- neutron star con- 
sisting of free baryons, electrons and positrons, and ac- 
tive neutrino-antineutrino pairs (and perhaps also muon- 
antimuon pairs and strange quark matter), all in ther- 
mal and chemical equilibrium at a temperature T ~ 
30 — 70 MeV and entropy per baryon s ~ 10. 

During the shock reheating epoch, at times post core 
bounce t p h ~ 0.06 — Is, neutrinos diffuse out of the core 



and deposit energy b ehin d the stalled shock, driving the 
supernova explosion [103]. 

In the r-process epoch neutrinos continue to diffuse 
out of the star at times i p b ~ 1 — 15 s. The neutrinos 
drive mass loss from the proto-neutron star, possibly set- 
ting the physical condition s for r-process (heavy element) 



nucleosynthesis 104 105 



The star's nuclear composition becomes increasingly 
neutron rich as the elec tron neutrinos depart the star 
and the reactions in Eq. ( 10.1 ) proceed with continuous, 
local re-establishment of (3 equilibrium. For example, at 
core bounce (t p b = 0), the electron fraction (net number 
of electrons per baryon) in the center of the core is Y e « 
0.35; a cold neutron star, the end point of evolution in 
this case, has Y e < 0.01. 

Since neutrinos dominate the energetics of core- 
collapse (Type Ib/c, II) supernovae, appreciable active- 
sterile neutrino transformation can completely alter the 
standard picture of stellar collapse. Although this picture 
has be en re fined by observations of neutrinos from SN 
1987A [106 1 , large-scale numerical simulations [103, 107 1 , 
and semi-analytical work [103, 108, 109, 110], it is not 



clear that this is the only way for supernovae to evolve. 
Nevertheless, by requiring that not too much energy 
be "lost" too quickly to singlet neutrinos in the proto- 
neutron star, and hence avoiding a conflict with stan- 
dard supernova theory and observations, we can delimit 
regions in the m s — sin 2 20 plane which may adversely 
affect its evolution. In fact, many supernova constraints 
on new physics rely on limiting the energy loss to exotic 
particles in a proto-neutron star 1 11 1 ] . 

Sterile neutrinos can be produced in supernovae either 
coherently through mass level-crossings or incoherently 
via scattering-induced wave function collapse. Whether 
the former or latter process dominates depends on the 
hierarchy of length scales relevant for neutrino trans- 
port and transformation: the local neutrino mean free 
path A, local neutrino oscillation length in matter l m , 
and resonance width Sr [the supernova analogue of Eq. 
(7.8)]. Neutrino flavor eigenstate evolution is coherent 
and transformation is uninterrupted through resonance 
if A ^> Sr. Neutrino flavor transformation is also adia- 
batic if Sr 3> /™ s , where Z™ s > l m is the oscillation length 
at resonance. This is in complete analogy to the evolu- 
tion through resonances discussed above in Sec. VIl for 
the early universe. If instead neutrino conversion is inco- 
herent and transformation is interrupted often by colli- 
sions, then the conversion will be suppressed by quantum 
damping if l m > A. 

It is easy to show that incoherent conversion domi- 
nates sterile neutrino production in supernova cores. In 
a proto-neutron star of radius R = 10 — 50 km the den- 
sity varies relatively slowly with distance from the center 
of the star, so the scale height of the weak potential is 
H = \d\nVfdr\~ 1 w 10 — 100km. Then the resonance 
width is Sr = H tan 20 » (10 6 -10 7 cm) tan 20. The mean 
free path of a typical neutrino is about A « 10 cm (actu- 
ally A ~ 10 cm to ~ 10 m, but this range makes little dif- 
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ference for our conclusions). Therefore, neutrino evolu- 
tion is coherent if sin 2 29 < 10~ 12 — 10~ 10 , and the maxi- 
mal resonance width giving coherent evolution is Sr max ps 
10 cm. Now the oscillation length at resonance is ZJ^ S = 
iirE/Sm 2 sin 29 ps 10" 4 cm(10keV/TO s ) 2 /sin26», so co- 
herent evolution gives Z™ s > 10 — 100 cm (10 keV/m s ) 2 . 
In these conditions, however, the evolution is not adia- 
batic, since <5r max < V^. As a result, mass level crossings 
contribute only subdominantly to sterile neutrino pro- 
duction in supernovae. Conversion to sterile neutrinos 
less massive than lOkeV is even less efficient, and more 
massive neutrinos will not have resonances. 

Incoherent neutrino production is weakly damped 
(or away from) a resonance when sin 2 29 



at 



> 



10~ 10 (10 keV/m 5 ) 4 . The local oscillation length in mat- 
ter away from resonance is typically much smaller than 
the oscillation length at resonance, so for sterile neutrino 
masses m s > 10 eV weak damping generally obtains off 
resonance even if the condition above is not met. Fur- 
thermore, since resonances last only a fraction < 10~ 8 of 
the neutrino diffusion time scale, as derived below, most 
neutrino conversion in a core-collapse supernova occurs 
in the weak damping regime. 

The time rate of energy "loss" £ to sterile neutri- 
nos (and sterile antineutrinos) per unit mass in a proto- 
neutron star is proportional to the active neutrino energy, 
scattering cross section on wea k ta rgets, and the average 
conversion probability 126 
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where, for the sake of simplicity, we have suppressed 
Pauli-blocking effects. In Eq. (10.2), the cross sections 
for neutrino scattering on free baryons are a Va (E) ps 
<jp a (E) ps l.QQG 2 F E 2 , and the differential neutrino and 
antineutrino fluxes are 
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(10.4) 



for relativistic neutrinos, where rj Va — fi Ua /T Ua is the 
denegeracy parameter of v a as above. A simple but 
somewhat crude criterion for avoiding conflict with su- 
pernova theory and observations of SN 1987A is £ <■ 
10 19 erg s _1 g _1 | 111 , a limit on the sterile neutrino emis- 
sivity equivalent to a loss of ~ lOMeV per baryon per 
second. 



The weak potentials [cf. Eqs. (5.5) and ( |5.6| )] which 
drive neutrino flavor transformation take a slightly dif- 
ferent form in supernovae than in the early universe flS4| . 
For v e v s transformation the p oten tial stemming from 
finite density effects is [from Eq. (5.5)] 
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(3F e - 1 + 4Y; c + 2Y V + 2Y Vt ) , (10.5) 



where m^r ps 931.5 MeV is an atomic mass unit, rij is the 
number density of species i, and = (n, — n-^jnB is 
the net fraction of species i per baryon (ng = n n + n p ). 



In part, the form of Eq. (10.5) follows from local charge 
neutrality and the assumption that muons, anti-muons, 
and strange quark matter are negligibly populated. The 



finite density 
[from Eq. (p. 
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transformation is 



V 



GfP 

V2', 



m N 



(Y e -1 + 2Y^ + 4Y I/(j + 2Y Vt 



(10.6) 



The poten tial for v T # v s transformation follows from 
Eq. ( |l0.6 ) upon switching the labels and v T . The 
finite temperature potential V T is negligible in super- 
novae (6^[, so the total potentials for neutrino and an- 
tineutrino transformation are V = V D + V T ps V D and 
V = —V D + V T ps —V D ~ —V, respectively. 

The average oscillation probabilities in Eq. ( 10. | ) de- 
pend on the weak pot e ntial s and neutrino mixing param- 
eters, from Eqs. ( p^ , (x5): 
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where A(p) = 6m 2 /2p ps m 2 /2E ps A(E) for rela- 
tivistic neutrinos. As before the quantum damping rate 
D = r„ Q /2 = J d^ Ua a Va {E)/2 for neutrinos is one-half 
the neutrino scattering rate; an analogous expression ap- 
plies for the antineutrino damping rate D. The coherent 
effects described in Ref. |69) may modify these conversion 
rates. 

The conditions for resonance ±V D = A (j p)cos 2(9 ps 
m 2 s cos29/2E, the potentials in Eqs. (|l0.5|), (|l0 6| ) and 



the average oscillation probabilities in Eqs. (10.7), ( |1 .S 
imply a negative feedback between active (anti) neutrino 
transformation and the potentials. If quantum damping 
effects arc comparable for neutrinos and antineutrinos or 
are relatively unimportant, as argued above for most of 
the parameter range of interest, neutrino conversion is en- 
hanced relative to antineutrino conversion when V D > 0. 
The preferential conversion of neutrinos decreases the po- 
tential, since the finite density part of the potential de- 
pends on the relative excess (or deficit) of particles over 
antiparticles (the net contribution of particle-antiparticlc 
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pairs is zero). As the potential decreases, the neutrino 
and antincutrino conversion rates approach the common 
rate they would have at vanishing potential. Since neu- 
trinos and antineutrinos of all flavors are produced in 
roughly equal numbers and their individual number den- 
sities are at least comparable to the electron density, the 
dynamical feedback mechanism may drive the potential 
V D locally to zero throughout the proto-neutron star, as 
long as the time scale for this process is less than the 
local neutrino diffusion time scale t d ~ X(R/X) 2 /c ~ 3 s; 
a similar sequence will ensue if V D < 0. If this happens, 
the equal neutrino and antineutrino conversion rates en- 
sure that the potential will remain zero, with any local 
deviation continuously smoothed out, up to the higher 
order effects of neutrino diffusion. 

In this sense zero potential is a fixed point of the 
dynamics of neutrino transport and conversion in core- 
collapse supernovae and the time scale for achieving zero 
potential is a dynamical equilibration time scale. The sit- 
uation is more complicated if sterile neutrinos mix with 
multiple active neutrinos or if damping is important. For 
example, an excess of neutrinos over antineutrinos im- 
plies damping is asymmetric, and this can retard equili- 
bration. These issu es will be examined in more detail in a 
separate work [ |112[ . We note in passing that this behav- 
ior can also occur when neutrino scattering is rare and 
mass leve l crossings dominate sterile neutrino produc- 
tion [113 1; some of this physics has also been discussed 
in connection with active-active neutrino transformation 



in proto-neutron stars [ 1 14 1 . 

We can estimate the equilibration time scale Ty for 
driving a potential v}f to zero by computing the time 
over which the disparity in the rates of neutrino and 
antineutrino conversion will induce a potential 5V D — 
— Vq° which cancels the original potential V® . Now Eqs. 
(goj), dial ) imply 5V D = 2y/2G F {8n Va - 6n Ba ), where 
8n Ua and Sn Pa are the numbers of neutrinos and an- 
tineutrinos converted per unit volume, respectively. The 
number densities of (anti)neutrinos converted are pro- 
portional to the rates of (anti)neutrino conversion per 
baryon, the baryon density, and the time over which the 
conversion takes place. If we assume Ty is small com- 
pared to the time scales for the change of the density and 
temperature in the star (a "static" approximation), we 



have 5n Ua = -n B Ty J d^ Ua a v ^(E)\(P m (v tl 



s-,P,t)} 



and similarly for 5n 9a . Solving for Ty and using rig — 
p/rriN gives the time scale for achieving dynamic equilib- 
rium in the v a ^ i/ s , v a ^ v s system: 



Ty 



V D m N 
2V2G 



1 



FP 



d^u a a Ua {E)-{P m {v a -> v s \p,t)) 



d$ Va <T Va {E)-{P m {D a ^v s -p,t))\ . (10.9) 



It is instructive to evalute the equilibration time scale 



Ty for the Vn ^ v. 



v s system in two limits: (i) far 



and (ii) very near or at a resonance, where 
A(£)cos20 w m 2 /2 E for 9« 1. In the o ff-res o nance 
case, inserting Eqs. fllO.ap-fllO.^) and Eqs. (|l0.7|)-([iolj|) 
in Eq. ( 10.9 ) and expanding to lowest order yields 



_off-r 



9G 



(VrP) 4 



m N 



lMV2G 3 F pT 2 m 6 s sin 2 29 
1.5 x 10" 8 s /10 14 gcm" 3 
sin 2 29 V P 
50 MeV\ 2 / 10 keV\ 6 / V D 



r ) 



leV 



(10.10) 



where the factor of 1.66 comes from the neutrino scat- 
tering cross section given earlier. We have derived Eq. 
( 10.10 ) by assuming a locally constant density and tem- 
perature, consistent with our "static" approximation. 
We have also ignored the build up of a non-zero muon 
neutrino chemical potential p v and accompanying Pauli 
blocking in neutrino scattering and pair production as 
equilibration proceeds, a simplification giving d$ U)i w 
d^ Ut . Using the same simplifications, we can find the 
time scale in the on-rcsonance case: 



on — res 



8tt 2 



\V D \m N 



(1.66)45^(5) G%pT 5 
16C(3) m 2 s m N 
(1.66)7tt 2 V2C(5) Gl-pT 6 

2xl0^s(^X) 5 ( 10l4 g Cm ' 3 



V D 



leV 



6.6 x 10~ 10 s 



50MeV 
f 



10 14 g cm 



-3 



lOkeV 



(10.11) 



The second and fourth expressions of Eq. (10.11) follow 
from the resonance condition « m 2 s /2E and taking 
the neutrino energy E to be the average neutrino energy 
(E) =fd$u ll E/fd§ Vlt » 7^/180C(3) w 3.15T„ M in 
the proto-neutron star. 

The on- and off-resonance equilibration time scales 
both vary inversely with the ambient density and tem- 
perature. Hotter conditions enhance the neutrino con- 
version rate per scatterer, and denser conditions en- 
hance the number density of scatterers. In either case 
it takes less time to drive a pre-existing potential to 
zero. On the other hand, only the off-resonance time 
scale T v s ~ res depends on the vacuum mixing angle; the 
inverse dependence is natural, since the emissivity in ster- 
ile (anti)neutrinos is proportional to sin 2 29. The on- 
resonance time scale r ° n_rcs is independent of sin 2 29, 
because the relevant mixing angle at a resonance is tt/4. 
If V® > neutrino conversion dominates antineutrino 
conversion at resonance and vice versa if VJj° < 0, so this 



away form a resonance, for which |V D | » A(E) cos 29; time scale is extremely short compared to the neutrino 
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diffusion time scaie or indeed any time scale typically as- 
sociated with a core-collapse supernova. Individual reso- 
nances are fleeting and result in negligible energy loss to 
sterile neutrinos. As a result, we may safely assume that 
most neutrino transformation takes place off resona nce, 
with an equilibration time scale given by Eq. ( 10.1 d| ) . 

We can delimit ranges of the neutrino mixing param- 
eters m s and sin 2 26 which may adversely affect core- 
collapse supernov ae by evaluating the sterile neutrino 
emissivity in Eq. (10.2) in the following manner. 

In the case of v e # v s transformation, the poten- 
tial V D is initially positive in the post-bounce supernova 
core. If the neutrino mixing parameters happen to give 
Ty ff-rcs > td, core neutronization proceeds roughly on a 
diffusion time scale, although the conversion of electron 
neutrinos accelerates the shift of (3 equilibrium toward 
neutron richness. Neutronization decreases the electron 
fraction Y e , and neutrino conversion decreases the elec- 
tron neutrino fraction Y Vb . As a consequence, the poten- 
tial decreases and eventually reaches a value such that 



off — res 



< td- 



'y ~- L d- 

Once this occurs, neutrino conversion quickly resets 
the potential to zero and maintains this value as neu- 
tronization continues. If instead the neutrino mixing pa- 
rameters happen to give Tv s ~ Tes < td at core bounce, 
conversion immediately resets and maintains the poten- 
tial at zero. 

In the case of # v s or v r ^ v s transformation, 
the potential V D is initially negative in the post-bounce 
core. Here conversion has no direct effect on the rate of 
neutronization (notwithstanding feedback on the nuclear 
equation of state) , which increases the magnitude of the 
potential as the electron fraction falls to a value ~ 0.01. 



Unless 



off — res 



< 



from the outset, neutronization dic- 



tates the evolution of the potential. 

The emissivity £ depends on the local, instantaneous 
value of the potential, so the total energy lost to sterile 
neutrinos depends on the time history of the spatially 
varying potential and, in particular, whether neutroniza- 
tion or neutrino conversion locally dominates its evolu- 
tion. An accurate estimate of the spatially integrated 
emissivity and, hence, of the regions in parameter space 
which may alter the standard picture of stellar collapse 
clearly requires a detailed investigation of the dynamics 
of these systems, as well as a better understanding of the 
physics of core-collapse supernovae. 

From the viewpoint of the viability of sterile neutri- 
nos as dark matter, however, we may conservatively esti- 
mate the emissivity by assuming the system achieves zero 
potential relatively quick ly, so that £ is evaluated with 
V D = 0. Evaluating Eq. ( 10.2 ) in this limit and imposing 
the condition £ < 10 19 ergs -1 g _1 gives the disfavored re- 
gions sin 2 29 > 3 x 10 -10 for u e v s transformation and 
sin 2 28 > 10~ 9 for ^ v s or v T ^ v s transformation for 
sterile neutrino masses m s > 10 eV. For smaller masses, 
the limits are significantly weaker, owing to the onset of 
quantum damping. These limits are shown in Figs. |^, 
p|, 0, and ^. The emissivity falls with increasing potential 



away from a resonance, so more detailed future work on 
this complex problem may find actual constraints which 
are weaker than the very conservative limits given here. 

Sterile neutrino dark matter parameters which lie near 
the edges of the disfavored regions in Figs. ||, ||, |L [?|, 
and H could give interesting signals in current and fu- 
ture supernova neutrino detectors (see, e.g., Ref. |L15| ). 
These detectors possibly could discern unique signatures 
for sterile neutrinos. Such signatures would bolster the 
case for sterile neutrino dark matter. 

To summarize this section, we have delimited conserva- 
tively the regions in the parameter space of active-sterile 
neutrino mixing that are disfavored by energy-loss con- 
siderations in core-collapse supernovae. We have found 
that the coupled problem of neutrino transport and fla- 
vor transformation in hot and dense nuclear matter is a 
formidable one. It involves following the local evolution 
of the weak potential which drives flavor transformation, 
including the feedback from diffusion and the conversion 
itself. We have described this physics roughly by esti- 
mating the competing time scales for lepton number dif- 
fusion and for cancellation of the potential. Depending 
on the mixing parameters and the spacetime evolution of 
the potential, the potential may well be reset to zero in 
much of the proto-neutron star, so the effects of neutrino 
propagation in matter need not suppress sterile neutrino 
(or antineutrino) production. As a result, core-collapse 
supernovae can be significantly more sensitive to active- 
sterile neutrino mixing than they were found to be in 
previous studies [^8|, [n], |(||. These studies used a spa- 
tially and temporally constant value of the potential, re- 
sulting in limits on mixing angles similar to ours but for 
significantly larger sterile neutrino masses > 10 keV; for 
smaller masses their limits on mixing angles are weaker 
because the putative matter effects suppress conversion 

Of course, we have deliberately chosen to be conser- 
vative in applying supernova limits, since our objective 
in this work is to assess the viability of sterile neutrino 
dark matter. The true limits, quite interesting in their 
own right and obtained from a self-consistent and proper 
treatment of the full, multi-dimensional Boltzmann evo- 
lution of the neutrino seas coupled with the nuclear equa- 
tion of state, may well lie somewhere in between the val- 
ues determined in this and previous studies. We leave 



attempts at such investigations for a future work [112 



XI. CONCLUSIONS 

We have estimated the resonant and non-resonant scat- 
tering prodution of sterile neutrinos in the early universe. 
The basis for the production of these sterile species is a 
presumed mixing with active neutrinos in vacuum. Of 
course, such mixing renders these species not truly "ster- 
ile." As a result, the "sterile" neutrinos can decay and 
this, together with their overall contribution to energy 
density, constitutes the basis for several stringent cos- 
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mological and astrophysical constraints which we have 
discussed in detail. 

Additionally, these sterile species may produce signif- 
icant effects in core-collapse supernovae. Some of these 
effects, such as massive core energy loss, could be the 
basis for true constraints. However, it must be kept in 
mind that (1) the supernova explosion energy is only 
some ~ 1% of the total energy resident in the active 
neutrino seas and that (2) we do not yet understand in 
detail how supernovae explode (nor do we have a suf- 
ficiently detailed observed core-collapse neutrino signal 
to place stringent constraints). Better theoretical un- 
derstanding of supernova physics, perhaps coupled with 
the neutrino signature of a Galactic core-collapse event, 
may allow the indicated regions on Figs. and || to be- 
come true hard and fast constraints instead of simply 
"disfavored" parameter regions. In fact, deeper insight 
into the time evolution of the potentials governing sterile 
neutrino production in the core may allow extension of 
the constrained parameter region to even smaller values 
of vacuum mixing angle. 

Nevertheless, it is clear from our work that sterile "neu- 
trino" species with ranges of masses and vacuum cou- 
plings could be produced in quantities sufficient to ex- 
plain all of the non-baryonic dark matter, while evading 
all present day laboratory and astrophysical constraints. 
Within the allowed ranges in mass and mixing parame- 
ters which give these dark matter solutions are regions 
where the sterile neutrino masses and/or energy spectra 
combine to produce collisionless damping scales corre- 
sponding to warm or even cold dark matter, subsuming 
the interesting behavior range for large scale structure. 

It is then a disturbing possibility that the dark mat- 
ter might not be "weakly interacting massive particles" 
(WIMPs), but rather "nearly non-interacting massive 
particles" (NNIMPs) which are likely not detectable in 
ordinary dark matter detection experiments. This pos- 
sibility begs two questions: (1) how could we hope to 
constrain or definitively detect this dark matter candi- 
date; and (2) what are sterile neutrinos? 

The answer to the first question is more straightfor- 
ward than the resolution of the second. As outlined 
above, better understanding of the neutrino and equation 
of state physics of core-collapse supernovae could help us 
extend constraints. Improved observations and models of 
x-ray emission from clusters of galaxies and other objects 
could provide stringent constraints or, conceivably, even 
detections of sterile neutrino photon decay. 

As discussed above, future observations bearing on the 
clustering of dwarf galaxy or Lyman-a cloud halos at high 
redshift may provide evidence for or against WDM as 
opposed to CDM. Direct evidence for WDM would con- 
stitute a point in favor of sterile neutrinos, though not a 
definitive one by any means. Likewise, and insidiously, 
our allowed parameter space for dark matter accommo- 
dates standard CDM behavior, even for lepton numbers 
equal to the baryon number. 

Direct detection of WIMPs in the laboratory or detec- 



tion of gamma rays asso ciated with WIMP annihilation 
in galactic centers [116] obviously rule out sterile neu- 
trino dark matter, given the small vacuum mixing angles 
suggested by our work. It is worth considering whether 
/3-decay electron energy spectrum or pion decay experi- 
ments could be pushed in sensitivity to the point where 
massive sterile neutrinos in some of the allowed regions 
of Figs. and || could be constrained. This would re- 
quire an increase in sensitivity to m s sin 2 29 of at least 
some six orders of magnitude and this is clearly untenable 
with current technology ]35| . Finally, although a number 
of extensions of the standard model motivate the exis- 
tence of multiple sterile neutrinos, it must be pointed 
out there is no independent physics suggestion for the 
sterile neutrino mass (~ I keV to ~ 10 MeV) and mix- 
ing (10~ 17 < sin 2 29 < 3 x 10~ 10 ) parameters which give 
viable dark matter candidates in our calculations. 

As discussed above, sterile neutrino degrees of freedom 
with ultra-large masses (e.g., of order the standard model 
unification scale, or even the top quark mass) are in some 
sense "natural," at least in the context of a see-saw ex- 
planation for the low masses of active neutrinos. 

There is now a reasonable chance that new neutrino 
experiments scheduled to come to fruition in the next 
few y ears [e.g., the Sudb ury Neutrin o Ob servat ory ( SNO) 



117) , KamLAND |ll8|, O RLaND JllJ, K2K [12C|, MI- 



NOS pl| , ICARUSp22|, and Mini-BooNE p3[ ] will 
allow us to deconvolve the neutrino mass and mixing 
spectrum. This would be an achievement heavy with im- 
plications for many aspects of physics and astrophysics. 
Will the requirement for sterile neutrinos remain? 

It is now true that the interpretation of the current so- 
lar, atmospheric, and accelerator (LSND) data in terms 
of neutrino flavor mixing physics demands the introduc- 
tion of a sterile neutrino with a rest mass comparable to 
that of some of the active neutrinos, i.e., light. Neces- 
sarily, then, this sterile neutrino is not the dark matter 
candidate we speculate on in this work. However, the 
unambiguous establishment of the existence of a light 
sterile neutrino would expose our ignorance of physics in 
the neutrino sector in a stark and dramatic way. On 
this score, the Mini-BooNE exp erim ent |123| and the 
SNO neutral-current experiment [117] are the most cru- 
cial ones for sterile neutrino dark matter. A confirmation 
of the LSND result invites speculation on the existence 
of more massive sterile neutrino states. 
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APPENDIX: THE GENERAL 
TIME-TEMPERATURE RELATION 

In this appendix, we review the calculation of the tem- 
perature evolution in the early universe through periods 
of varying statistical weight in relativistic particles, g, 
Eq. ([7^). The time derivative of the temperature can be 
written as [ 1 24 1 




(1) 



where r = ln(i? 3 ) and R is the scale factor. The expan- 
sion rate is determined by the Friedmann equation 



H=§ ^ = — V^ot- 0.207 s-V/ 2t2 , (2) 
at K m p i V o 

where p to t = (tt 2 /30)gT 4 is the total energy density, and 
T is in MeV in the last approximation of Eq. (||). There- 
fore, the first half of the time-temperature relation, Eq. 
(|l|), is straightforward: dr/dt — 3H. One can get the sec- 
ond half through the conservation of comoving energy, 



d 



dt 



This can be rewritten into the desired form 



dr 
dT 



rfptol 

dT 



{ptot + Ptot) 1 



(3) 



(4) 



With Eqs. @ and (§) one has the general temperature 
evolution. 

In our case the sterile neutrinos can contribute sig- 
nificantly to the energy density and pressure. Approxi- 
mating all species other than the sterile neutrino to be 
relativistic, we havep* « 1/3/9*, where and are the 
energy and pressure in all particles other that the sterile 
neutrinos / antineutrinos. Therefore, 



dT 
dr 



dp* 



dp s 
dT 



(5) 



The rate of change of the standard energy density is 
straightforward: 



dT 30 9 30 dT' 



(6) 



The temperature derivative of the sterile neutrino energy 
density is 



dp s 
dT 



{U, (P) + A (P)) [{P 2 + m 2 s ) /T 2 ] 1/2 p 2 dp 



$ J IfvM + fM) [{p 2 +m 2 s )/T 2 ]- 1/2 p 2 d P . 



This, together with the energy density and pressure in 
the sterile neutrinos and antineutrinos calculated from 
their time-dependent distribution functions, allows one 
to readily arrive at a consistent time-temperature evolu- 
tion. 
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